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The interaction between planetary waves and an arbitrary zonal flow is studied using wave-kinetic 
methods. A quantum mechanical tool, the Wigner distribution, is used to construct a phase space 
transport equation for the waves that does not depend on any assumption of scale separation, and it 
is shown that this transport equation reduces to a standard wave-kinetic form when scale separation 
can be assumed. Analytical arguments, combined with a simple numerical model, are then used to 
demonstrate that at least in the weakly turbulent case, both zonal jet formation and east-west jet 
asymmetry can be explained intuitively in terms of the motion of planetary wavepackets in phase 
space. 



I. INTRODUCTION 



The formation and maintenance of large-scale jets by the collective nonlinear interaction of eddies is a fundamental 
problem in planetary fluid dynamics. Perhaps the most famous example of the effect can be seen in the atmospheres 
of the planets Jupiter and Saturn; it has long been believed that the stable, coherent jet structures observed there 
owe their existence to continual forcing by smaller-scale eddies. Other geophysical examples of interest include the 
Earth's atmospheric jet stream and, quite possibly, the alternating zonal jets recently observed in the Earth's Pacific 
ocean [§]. 

In all these cases, the change of planetary vorticity with latitude or /3-effect is believed to be a vital part of the 
problem, as it allows for the presence of planetary waves in the system. In a seminal paper, Rhines studied the 
interaction between planetary waves and turbulence, and came to the conclusion that the inverse energy cascade of 
'ideal' two-dimensional turbulence would be halted by planetary wave motion at large scales, leading to the transfer 
of energy into the zonal modes and hence to jet formation. Although extremely insightful, Rhines' work was partly 
heuristic and could not provide a detailed dynamical explanation of the process. 

As a result, the theoretical investigation of reduced problems involving the interaction of zonal jets and planetary 
waves is still of great importance to our overall understanding of atmospheric and oceanic fluid dynamics. Much 
interesting work has previously been done on the subject; for example, a number of authors (|13j. [16] . [6]) have 
studied the interactions of planetary waves with critical layers in various idealised scenarios. The foundation for 
many of these studies was the earlier development of various real space conservation theorems (see e.g. [1], [5]), the 
most important of which are now regarded as an essential part of basic wave - mean flow theory. 

Other researchers, beginning with Dyachenko et. al. [5], have taken an approach based on wave-kinetic methods. 
Among these were Benilov et. al. [3J, who studied the interaction of planetary waves with a single jet and noted the 
inhomogenous, antifrictional properties of the system, and Manin and Nazarenko |5], who used a wave-kinetic equation 
to study the interaction between scale-separated zonal flows and planetary waves in the limit of small /3-effect. Most 
recently, there have also been some attempts to utilise the wave-kinetic approximation in wave - mean flow numerical 
simulations (see e.g. [7] and [T5]). 

In this work the aim is to combine the two approaches. For the first time, techniques are developed for the study 
of the interaction between an arbitrary, finite amplitude zonal flow and a broadband distribution of planetary waves. 
The Wigner distribution, a tool first developed for use in quantum mechanics |17j . is used to describe the wave-mean 
flow interaction purely in terms of wavepacket motion in phase space. As an example of the applications of these 
techniques, we then study the onset of jet formation and asymmetry due to the nonlinear interaction of a planetary 
wavepacket with the zonal flow. 

In Section [IT] existing wave-mean flow theory for quasigeostrophic flows is briefly discussed. In Section |III| the 
Wigner distribution is introduced and used to derive the 'quantum' Boltzmann equation that describes the transport 
of planetary wavepackets, and it is shown that an assumption of scale separation leads to a wave-kinetic equation of 
the type discussed by [5 and others. Finally, in Section IV a simple numerical simulation is introduced and is used 
to demonstrate the utility of the Wigner distribution in analysing wave-mean flow problems. 
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FIG. 1: Schematic of the theoretical setup: a /3-plane model periodic in the z-direction and open in the {/-direction. The /3-plane 
approximates fluid motion on the midlatitudes of a planet, with x and y equivalent to east-west and north-south directions 
respectively. 



II. FUNDAMENTALS OF WAVE-MEAN FLOW THEORY 

The essential features of many large-scale geophysical flows can be captured by the quasigeostrophic potential 
vorticity (QGPV) equation 

= § + = (D 

where ip is the velocity streamfunction, q = (d xx + d yy + d z ((fo/N 2 ) d z )) i\) + py is the quasigeostrophic potential 
vorticity (3 is the local gradient of planetary vorticity with latitude, J[, ] is the two-dimensional Jacobian operator 
such that J[f, g] = d x fd y g — d x gd y f in Cartesian co-ordinates, and k is the Ekman damping parameter. Equation 
is simply a statement that q is conserved following fluid elements in the absence of damping and forcing. It is an 
approximation to the full Navier-Stokes equations that is applicable when the system under consideration is rapidly 
rotating and strongly stratified, and is derived in detail in many good fluid dynamics textbooks (for example see |12j). 
For convenience, in this paper we will work with ([I]) in Cartesian co-ordinates according to the standard /3-plane 
model (see Figure [lj . We also focus on the situation where the system is unbounded in the y (north-south) direction. 

As a result of the background vorticity gradient [3, the linearised form of has planetary wave solutions, with 
dispersion relation 
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when no zonal flow or damping is present. In ([2]), k z is understood to be an eigenvalue of the usual vertical structure 
equation such that d z ((fo/N 2 ) <9 2 <I') = —k 2 ^, subject to suitable boundary conditions. 

To investigate the interaction between planetary waves and zonal flow, it is standard to define an average in the 
x-direction such that any quantity decomposes into a mean flow and disturbance field: /(x, y, t) — > f(y, t) + f'(x, y, t). 
Then, ([!]) becomes 

m=-dy Vq Kq (3) 



for the mean flow and 



dq' 
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for the disturbances, with 7 = d y q the total gradient of zonal potential vorticity. The left hand side of Q describes 
the evolution of planetary waves in the presence of a zonal flow, while the right describes nonconservative effects and 
wave-wave interactions. 

As this work is primarily concerned with wave-mean flow interaction, we will assume the wave- wave terms in Q to 
be small from here onwards. In Section [TV] where we discuss the numerical simulation of a generic jet- wave interaction 
problem, the situations where this assumption begins to fail will be made clear. 

For moderate disturbance amplitudes, it can be shown that the wave action, defined as be \q' 2 /l, is a conserved 
quantity. This can be seen through multiplication of Q by q' which results in 

Qji 

— + v'q' = -2/cn. (5) 
at 
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if terms of order q' 3 and greater are neglected. For further details of the derivation, see pQ. If scale separation between 
zonal flow and waves is then assumed and the waves are taken to be monochromatic, a real space transport equation 
for wave action can be written 

-Qj. + V 2 d ■ (v 2 £>n) = -2/cn, (6) 

where V2D = (d y ,d z ) and V2D is the (y,z) group velocity of the waves (see for example [ID] , page 403). 

The final established result of importance to the rest of this paper is the powerful non-acceleration theorem, which 
states that in the absence of forcing or damping, the rate of change of zonal potential vorticity and wave action are 
directly tied to each other 

£(*-*)-* n 

Equation ^ can be proved by use of the Taylor identity [I], or by combining equations ([3| and ^ and setting k — 0. 



III. THE 'QUANTUM' BOLTZMANN EQUATION 

Having described the basic features of the planetary wave-mean flow problem, we now wish to extend the scope of 
the analysis. The aim is to construct a generalised phase space transport equation for the planetary waves, leaving 
any scale separation approximations until as late as possible in the derivation. 

Perhaps the best known example of a phase space transport equation is the Boltzmann equation of classical me- 
chanics. It describes the (smoothed) collective interaction of an ensemble of point particles in phase space, due to a 
field that is a function of the position of every particle. For example, in stellar dynamics, the Boltzmann equation is 
used to describe the motion of an ensemble of stars, with each star's motion determined by the mean gravitational 
field due to all the others. 

In quantum mechanics, the quantum Boltzmann equation is needed to describe the collective interaction of an 
ensemble of particles. Due to the uncertainty principle, quantum particles cannot simultaneously have a well-defined 
position and momentum, and hence do not occupy single points in phase space. 

As fluid mechanical wavepackets are not pointlikc particles in phase space cither, if we wish to construct a general 
equation for their motion, we must initially take the quantum-mechanical route. The central difference between a 
quantum system and a fluid one is that in the fluid case, any wavefield can be decomposed into any number of 
wavepackets, as the base packet amplitude can be chosen arbitrarily. While wave action, a continuous quantity 
analogous to wavepacket number density, is conserved, the discrete notion of 'number of wavepackets' is nonsensical, 
strictly speaking, in a classical fluid context. 

In the case of quasilincar planetary waves, interaction can only occur indirectly through modification of the zonal 
flow profile. It is appropriate for this problem, therefore, to view the zonal flow as the 'mean field' through which the 
waves couple, analogous to the gravitational field in the stellar case. 

We begin the derivation of the planetary wave Boltzmann equation by writing the disturbance equation Q in terms 
of a new variable <j> = q' /y/7rj. If we neglect terms of order 2 and higher, Q then takes the form 

<%-**• < 8 > 

where the wave operator G is defined as 

_ —k _ 

G[x, k, t] = V7 jz ~* - f Ji + uk x -in (9) 
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and the position and wavevector operators are 

x = x and k = (-id*, -id y , -i^/dJJJi/N^ydSj (10) 

respectively, as we are working in a position space representation. Note that the denominator in Q is simply the 
potential vorticity inversion operator, such that 

- (kl + k 2 y + kfj' 1 q = (d XX + dyy + ^(/oViV 2 )^)- 1 q = V- (11) 
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FIG. 2: A real space function 4>{y) (left), its normalised Fourier power spectrum (bottom right) and the magnitude of 

the associated Wigner distribution (top right). In this example lo = 50, yo = 0.5 and Ay = 0.2. 



When scale separation is assumed, operators become real numbers, and the large-scale zonal flow 'sees' wavepackets 
as phase space points with exact values of x and k. Then, ^ simply becomes the generalised dispersion relation for 
small-scale planetary waves in the presence of a zonal flow and damping 

G[x, k, t] -> uj(y, z, k, t) = ^ x + uk x - in. (12) 

K x + K y + K z 

To progress towards the wave-kinetic equation, we now need to utilise a tool from quantum mechanics: the Wigner 
distribution. It is defined in three dimensions as 

A/^(x,k,t) = -^3 J "^(x-ixi^e-^^x+ixi^dSxj (13) 
or alternatively in spectral space as 

l f + °° (. l. A ^ . (. l 
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AT^(x,k > *) = ^3 / 0>* (k--k 1 ,i)e lxkl 4>(k+ik 1 ,t)d 3 k 1 (14) 



where $(k,t) = (2tt)~ 3/2 /+" exp[-«k • x]<£(x, t)d 3 x is the Fourier transform of <f>. A/^^ can broadly be thought of 
as a phase space distribution for the function <j>, but it has some fairly weird properties — not least of which being 
that it can take negative values. However, its projections onto real and spectral space are always positive valued. 
As an example, let us examine the simple wavepacket function 

(f>(y) = e -((y-yo)/Ay) 2 +U y ^ 

In this case, the Wigner function only varies in (y, k y ) space, where it takes the form 

A/^(y, k y ) = i J + °° <j>* (y ~ l -y^j e~^4> (v + \v^j *Vi (16) 
= (Ay/\/2~7r) e - 2 ((y-yo)/A a ) 2 -((fe H -i )Ay) 2 /2 ^ 

i.e. that of a two dimensional gaussian. Both <f>(y) and Af^^iy, k y ) for this example are plotted in Figure [2j along 
with the Fourier transform amplitude 

By taking a time derivative of (13) and using ([IJ, we may write 

( 18 ) 

Then, by defining the phase space operators X = x + | Vk and K = k — |V X (see e.g. [14| ) and noting that 

= -A/^x* = A/^, (19) 
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and hence clearly 



(20) 



the fairly weak assumption that G[X, K, f] can be expanded in powers of the two operators X and K allows us to 
arrive at 



i~=(G[X.K.l] -(V[X ! .K'./]).V. 



(21) 



where from now, Af = Afs,6 is written for brevity. 

Equation (21 1 is the 'quantum' boltzmann equation, describing the transport of planetary wave action in phase 
space. Although it does not depend on any assumption of scale separation between zonal flow and waves, we can 
easily show that it reduces to a more familiar wave-kinetic equation via a Taylor expansion of the operator G about 
x and k. Truncation of the Taylor expansion at first order allows us to write 



G[X,K,t] Ma;(x,k,t) 



dui i d 
dx ' 2dk 



dui i d 
5k ' 2dx' 



A similar expansion for G[X*,K*,t] and substitution into (21 1 then leads to 

I - V- =T[N] 



dM 



dM 
<9x 



dM 
9k 



(22) 



(23) 



where v and F, the group velocity of, and force on, a wavepacket respectively, have their usual definitions as v = Vk^ 
and F = — V x cj. Equation (23 1 is equivalent to (211 in the geometrical optics limit of small-scale disturbances. It 
describes the collective motion of an ensemble of point- like wavepackets through phase space. The right hand side of 
(23 1 contains all nonconservative terms: according to our derivation, T[M] = —2kM. However, if the effects of wave- 
wave interactions were to be included, T would also contain more complicated terms describing collisions between 
wavepackets (see for example [S])- 

For later use, we note that for small-scale planetary waves, we can differentiate (12 1 to obtain group velocities and 
rates of change of wavenumber or 'forces' 



rfx 
~dt 
dk 
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(24) 



where we have written |k| 2 = fc 2 , + ky + k 2 z etc. Note finally that if we define the wave action n (see equation (j5|) to 
be the projection of Af onto (y, z) real space 



/+00 
AA(x,k,t)d 3 k, 
-OO 



(25) 



we can recover the results of standard wave-mean flow theory outlined in Section [TT] By integrating ( 23 ) over k, 
making use of the fact that V x • v + Vk • F = and assuming that M — > as |k| — > 00, we arrive at 



dn 

at 



+ V 2 d • ((v 2 d)«) = -2/sn, 



(26) 



where 



+ OO 



vAAd 3 k 



(27) 



is the phase space averaged (y, z) group velocity for the planetary waves. Equation (26) is a generalisation of ^ to a 
broadband distribution of small-scale waves. It is interesting to note that its derivation from (23 1 is directly analogous 
to the derivation of the equation of continuity from the Boltzmann equation in fundamental fluid dynamics. 
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IV. NUMERICAL SIMULATION 



We now wish to put the Wigner distribution to work by considering a simple wave-mean flow numerical simulation. 
The essential features of the problem are captured by restricting the planetary wavefield to a single x-wavenumber, 



k x — ko, but allowing it to be broadband in k y . This is justified by noting that according to (24 1, the zonal mean flow 
cannot move the planetary wavepackets to different k x and also that in a wave-mean flow context, their x-position is 
clearly irrelevant. For simplicity we also restrict the problem to the single layer barotropic case, although it should 
be noted that all of the theory presented in the previous section is also applicable to mixed barotropic / baroclinic 
flows, which in general will vary with height as well as latitude and longitude. 

For barotropic planetary waves, only the y and k y dimensions of phase space are of relevance. In particular, the 



equation for phase space velocity vectors ( 24 1 simplifies to the two components 



_ 2jk X ky F ^j^ £> l_ k ^ (OK) 

Vy ~ |k|4 ' ry - |k|2^ x d y - ( ' 

To investigate equations ^ and (jJJ) numerically, the disturbance vorticity was written as q' = Ke[Qe th " x ], allowing 
the derivation of the simplified equations 

dO / \ d 2 

i^ = k (uQ + 1 L- 1 Q)-i K Q L = __fcg (29) 



and 



^ = v'q* -ku=~y (lm[i- 1 g]Re[Q] - Im[Q]Re[i _1 Q]) - ku (30) 



for waves and mean flow respectively. For all the numerical results presented here, equations (|29j) and (30 1 were solved 
using an explicit 4 th order Runge-Kutta method. The code used was designed to halt whenever a) the Rayleigh-Kuo 
criterion f3 — u" < for barotropic instability or b) the wave-breaking criterion |t/| TOO:r > u>/k x were satisfied. This 
ensured that the original physical assumptions behind the model were not broken during the simulation. 

As the simulation described is highly idealised and not intended to directly model real planetary flows, dimensionless 
units are used throughout this section. For comparison, however, we note that for a midlatitude slice of Jupiter's 
atmosphere, when scaled into units of planetary rotation period Tj and radius rj, the mean zonal wind speed is 
approximately u — 0.01 rj Tj _1 and the (3 parameter is /3 = 5 — 10 rj _1 Tj , depending on latitude. For 
all simulations here we used j3 = 10, and maximum zonal wind speeds of barotropically stable jets were of order 
max [it] = 0.001. Thus we are investigating a fluid dynamical regime with slightly weaker zonation, generally, than 
that observed on the gas giant planets. 

Firstly, we study the extremely simple case of an near-infinitesimal wavepacket with no initial zonal flow and no 
Ekman damping. Initial disturbance vorticity is 

Q = Q exp [il y - (y - y ) 2 /(Ay) 2 ] , (31) 

with Iq = k — 60, ?/o = 0.25 and Ay = 0.1. In Figure [31 the magnitude of the phase space distribution l-A/^l is 
plotted above zonal velocity u for a series of timesteps. As can be seen, when the wavepacket has wavevector such 
that k x k y > 0, it drifts northwards due to the /3-effect. Weak zonal jets form as a result of this motion. By the 
barotropic version of the nonacceleration theorem ([7| 

— - -— (32) 
dt ~ dt {6l) 

we see that latitudinal planetary wavepacket motion must always cause jets to form in this way. The process is 
demonstrated schematically in Figure [4] 

The group velocity calculated from p8| agrees closely with the observed velocity of the wavepacket peak (the 
difference is less than 3% in the example shown). However, note the stretching of the wavepacket in phase space. This 
is due to the dispersive nature of the planetary waves, as determined by Q. Essentially, the group velocity v v on the 
left hand side (in phase space) of the wavepacket is greater than that on the right — this is shown schematically by 
the arrows on the first plot in Figure [3] 

The second basic case of interest is that of an infinitesimal wavepacket on a linearly sheared zonal flow of the form 
u = — A(y — j/o)- Here, A = 0.01, yo = 0.5 and all other parameters are as in the previous example. As shown in Figure 
[5] in this situation a wavepacket with k x > is forced towards higher k y wavenumbers, losing energy to the zonal 
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FIG. 3: A planetary wavepacket on a /3-plane with positive k y and k y will move northwards. As group velocity depends on 
k y , the wavepacket becomes tilted in phase space (top row), although its volume remains approximately the same. Note the 
small-amplitude zonal flow (bottom row) induced by the wavepacket motion. 




FIG. 4: Schematic explanation of the jet formation seen in Figure [3} If a planetary wavepacket is moving northwards such that 
in time St = t% — t%, Sn — n(t2) — n{t\), then the nonacceleration theorem (32| ensures the zonal flow produced 8u will be of 
the form shown. 



flow in the process. As = 7 in this example, the enstrophy of the wavepacket remains constant and hence energy is 
transferred upscale, while enstrophy moves downscale. Again, scale separated predictions match the observed value 
closely for this case. 

More interesting and subtle phenomena occur when we allow the wavepacket to be of large enough amplitude for 
coupling with the zonal flow, but not wave-breaking, to occur (see Figure [6]). Then, we expect it to initially move 
northwards, with two zonal jets forming due to the non-acceleration theorem as in the first example. However, as the 
zonal flow becomes stronger, it begins to influence wavepacket propagation through a) the shear effect described in 
the second example and b) alteration of the basic potential vorticity gradient 7 = f3 — d yy u (see equation ( 28 > ) . As the 
initial zonal flow gradient in the centre of the channel is negative, the wavepacket is forced to higher k y wavenumbcrs, 
reducing its group velocity and hence the growth rate of the zonal flow. This process continues until the zonal flow 
either removes most of the wavepacket energy and reaches a quasi-steady state, or sharpens to the extent that it 
becomes barotropically unstable. 

It is fascinating to note that east-west jet asymmetry occurs regardless of the wavevector k sign of the initial 
wavepacket. As expected from (28 1 and shown in Figure [7j if the product k x k y is negative, the packet initially moves 
south, and the initial induced jets are of opposite sign. However, whatever the sign of k x and k y , jet formation always 
pushes some wave action to higher absolute wavenumber values \k y \. Combined with the fact that the wavepacket 
propogates away from the jet in a direction dependent on k x k y , the result is that in each case, the eastward jet 
becomes sharper than the westward one. 

This phenomenon can also be explained via a conservation of energy argument; jet formation increases zonal flow 
energy, which necessitates energy loss by the wavepacket. As the zonal flow cannot alter k Xl this can only be achieved 
through transport of wave action to higher absolute k y . Thus, east- west asymmetry appears to be a quite generic 
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FIG. 5: In the presence of a shear flow that decreases linearly with latitude (far left), a planetary wavepacket with positive k x 
will move toward higher wavenumbers, losing energy to the zonal flow in the process. Note the slight northward drift of the 
wavepacket; this is due to the /3-effect shown in Figure [3] 



time = time = 50 time = 100 time = 150 
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FIG. 6: Nonlinear evolution of a wavepacket (Qo = 0.12, j/o = 0.5, Ay = 0.05, all other parameters as in first example). 
Initially, the wavepacket moves as in Figure [3] but the zonal shear it produces modifies its motion as time progresses. Note 
that by t = 100, the Wigner distribution has become negative-valued in places. 



feature of wave-forced jets on the /3-plane that does not require the existence of critical layers for its occurrence. The 
fundamental cause of the phenomenon is evidently the chiral (left / right hand) symmetry breaking due to the fact 
that we are working in a rotating coordinate frame. This is manifested in the quasigeostrophic equation ([lj through 
the positive- valued (3 parameter. 

Use of the Wigner distribution allows us to examine in detail the regions where the wave-kinetic approximation 
(equation (23)) does badly, at any time during the simulation. In Figure [8j for an example nonlinear wavepacket in 
phase space(initial conditions Iq = = 30, yo = 0.5, Ay = 0.05 and Qo — 0.3), the actual time rate of change of 
the Wigner distribution dN^^/dt at t = 50 has been plotted alongside the wave-kinetic estimate of rate of change, 
which was calculated numerically for a single timestep according to ( 23 1 . The agreement between the two quantities is 
good near the centre of the wavepacket, but serious deviations begin to occur in the lower wavenumber region, where 
rapid fluctuations in the distribution are occuring (bottom left of the centre and right hand side plots) . Interestingly, 
however, for many cases where wavepacket - zonal flow scale separation could not be assumed, the wave-kinetic 
estimate of dAf^/dt continued to agree surprisingly well with the real value. 

If Ekman damping is present and the planetary wavepacket is continually forced, in many cases the system will 
eventually settle into a steady state. Figure [9] shows a plot of wavepacket and zonal jet total kinetic energy as a 
function of time for a simulation with k = 0.01 and an extra term added to the right hand side of (29 1 of the form 
+inQf, such that the disturbance quantity Q is continually relaxed to the profile Qf. For the simulation shown, 
Qf = 1 x 10 -1 / 2 exp[il y — (y — y a ) 2 /(Ay) 2 ]. It can be seen that after a period of rapid initial wavepacket growth, the 
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FIG. 8: Magnitude of Wigner distribution lA/^.^l (left) plotted alongside actual local time rate of change dAT^^/dt (center) 
and scale-separated dN /dt (right) for a generic wavepacket motion problem. The same contour scale has been used for the 
second two plots. 



zonal flow begins to gain energy by moving wave action to higher wavenumbers. This causes the total eddy energy to 
decrease, and hence the total system energy then slowly decreases until a steady state is finally reached. 
For a small-scale wavepacket, we expect the steady state zonal velocity to be given by 



v'q' 



d 
dy 



4-oo 



M{y,ky)vydk y 



(33) 



where v y is the group velocity in the y-direction and J\f(y,k y ), as defined in (16 1, is the projection of the Wigner 
distribution (13) onto two-dimensional (y,k y ) phase space. 

Now if the amplitude of the steady state wavepacket is small, we expect it to be relatively undistorted by the zonal 
jets it produces, and hence approximate the standard form 



Af(y, k y ) = ^ Qe -^(y-yo)/^y) 2 -((k y -io)A y r-/2 



(34) 



where Ao is the steady state peak wavepacket value. Provided that the variation of group velocity with k y is small 
over the width of the wavepacket, ( 33 ) can then be integrated to give an estimate of the total kinetic energy of the jet 
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FIG. 9: Eddy, zonal and total kinetic energies as a function of time for a simulation with planetary wavepacket forcing and 
Ekman damping of both zonal flow and eddies. 




FIG. 10: Total steady-state jet kinetic energy as a function of scale separation for the relaxation problem described in the text. 
Solid and dotted lines show linear analytical predictions and numerical results, respectively. Eddy forcing amplitudes were 
Qo = 1 x 10 -3 / 2 , 1 x 10" 1 and 1 x 10~ 1//2 for the left, centre and right hand side plots respectively. The right hand side plot 
has fewer data points because wave-breaking begins to occur at high forcing amplitudes. 



We expect (351 to be valid in the limit of forcing by small amplitude, scale-separated wavepackets. In Figure 



10 we have plotted the predictions of (351 against numerical jet kinetic energy for varying values of wavepacket 



scale seperation, denned as lo/k&y. Calculations were made for three different values of forcing amplitude, and each 
individual simulation was run for 100 times the relaxation time n~ l to ensure that a steady state was reached. For 
the plots shown, the east-west wavenumber is fixed at a constant value k = 40. As can be seen, the final jet kinetic 
energy rapidly increases to a maximum at around lo/k&y = 1, and then decreases, as the forcing wavenumber is 
increased. 

Note that as expected, the numerical results closely match the linear predictions only when wavelength is small 
compared to the jet width. In general, the analytical predictions tend to overestimate total zonal kinetic energy at 
low values of of lo/k& y . In the rightmost plot of Figure 10 the forcing amplitude is great enough for wave breaking 



to occur in some cases; this is the reason for the missing aata points in the centre of the plot 



V. CONCLUSIONS 



In this paper, the nonlinear interaction between an arbitrary zonal flow and a broadband distribution of planetary 
waves has been studied. The Wigner distribution has been used for the first time to derive a phase space planetary 
wave transport equation, (21), that does not depend on any assumption of scale separation. It has then been shown 



via a Taylor expansion that (21 1 takes the form of a standard wave-kinetic equation in the geometric optics limit. 

A simple numerical model has been used to investigate the onset of jet formation in a generic zonal mean - planetary 
wave interaction problem. By studying the motion of the planetary wave action in phase space via the the Wigner 
distribution, it has been possible to intuitively explain the eventual asymmetry of the fully developed jets. Also, the 
steady-state properties of the system in the presence of forcing and Ekman damping have been studied and it has 
been found that while linear analytical predictions hold well for small zonal jet amplitudes, they overestimate total 
jet energy when nonlinear effects become important. 

It is interesting that wave-kinetic arguments, such as those based on the ray-tracing equations (28 1), sometimes 
continued to hold fairly well even when scale separation could not strictly be assumed in the simulation. Possible 
extensions of ray-tracing ideas to more general systems have been noted before (e.g. |T|, page 188), and it is possible 
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that the subject deserves more thorough investigation. The work here opens up a way to do this analytically, for 
example by progressively including higher order terms in the Taylor expansion of equation (21 1. The problem could 
also perhaps be tackled numerically, by comparing a general wave-mean flow simulation of the type described in 
Section IV with a scale-separated one based on equation ( 23 1 . 

The main limitation of the analysis in this paper has been the restriction to wave-mean flow interactions only. 
Wave-wave and turbulent interactions can play a very important role in the overall development of fluid flows on 
the /3-plane, and any eventual general theory should aim to take them into account. To this end, it would be most 
interesting to extend phase space analysis techniques to a study of more strongly nonlinear planetary-scale turbulence. 
As a first step, such a study would require the derivation of a wave-kinetic equation for fully coupled planetary waves 
in the presence of an arbitrary zonal flow. As discussed by [8] and others, there are difficulties involved in this, as 
a standard derivation requires knowledge of the planetary wave normal modes, which for the general disturbance 
equation (4| have never been derived analytically. However, it is possible that progress could be made by taking an 
approach based on the real space definition of the Wigner distribution (13 1. 

More generally, it appears standard in nearly all studies of turbulent flows to take cither a real space view of 
the dynamics, or a spectral one. Both approaches have their respective advantages and drawbacks, but phase space 
analyses, by use of the Wigner distribution, allow us to take both simultaneously. Hence it is suggested that for many 
fluid dynamical problems, particularly those involving a mean flow and a disturbance, the techniques described here 
have the potential to greatly increase our understanding of the system's behaviour. 
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